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ABSTRACT 


In this thesis, a data fusion problem when a number of different types of sensors 
are deployed in the vicinity of a ballistic missile launch is studied. An objective of this 
thesis is to calculate a scoring function for each sensor track, and the track file with the 
best (optimum) track score can then be used for guiding an interceptor to the threat within 
the boost phase. Seven active ground-based radars, two space-based passive infrared 
sensors and two active light detection and ranging (LIDAR) sensors are used to track the 
ballistic missile in the boost phase. Each space-based platform carries one passive 
infrared sensor and one LIDAR. For the threat scenario, an IMPULSE intercontinental 
ballistic missile model is used to create the trajectory of a generic ballistic threat. The 
IMPULSE model is developed by the National Air and Space Intelligence Center to 
provide an accurate representation of ballistic missiles. Each sensor provides a track of 
the missile in the boost phase by using a multiple hypotheses tracking algorithm with an 
extended Kalman filter. The calculation of the track scoring function is to identify the 
sensor with the best track file. A track score is calculated for each sensor based on the 
kinematics of the missile flight parameters and the signal-to-noise ratio at the sensor. By 
using likelihood ratios, the optimum track file of the threat can then be determined and 
the corresponding track file can be transmitted to the battle manager control in order to 
lead the interceptor vehicle against the threat using the track file with the best score. 
Using the optimum track file scoring signal processing techniques developed in this 
thesis, the best track file can be sent to the interceptor to destroy the ballistic threat. This 
leads to a faster response management where the threat can be destroyed inside the 
territory of the country which launched the threat before any countermeasures are 


deployed. 
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EXECUTIVE SUMMARY 


During the past twenty years — specifically after the first war in the Persian Gulf 
where ballistic missiles were used for first time after World War I] — the global 
community started to realize the danger of ballistic missiles. An example is the case 
where these weapons are controlled by terrorist groups or imperialistic governments. 
Western governments, having realized the magnitude of the problem, started to develop 
programs with different types of sensors placed on different platforms in order to create a 


shield against the ballistic missile threat for their territories and their allies. 


The main problem for the deployment of all these sensors under the control of a 
central battle manager control is the evaluation of the sensor’s results and the subsequent 
analysis. The battle manager control needs to know in only a few seconds the exact 
position of the threat missile and its profile in order to send the appropriate interceptor. 


The goal is the destruction of the missile within the boost-phase of the missile flight. 


The geographical area of the scenario chosen for this thesis’ engineering analysis 
and program development is chosen without any political motivation. Also the type of the 
threat missile that is used in the study is an unclassified experimental model that is 


generic. 


For the threat scenario, an IMPULSE intercontinental ballistic missile (ICBM) is 
used to generate an attack against a distant country. The IMPULSE model is developed 
by the National Air and Space Intelligence Center to provide an accurate trajectory of 


ballistic missiles. 


In this thesis seven active ground based radars are deployed in the vicinity of a 
ICBM launch site in order to monitor its activity. Also, two low earth orbit (LEO) 
satellites containing passive infrared sensors (IR) and active light detection and ranging 
(LIDAR) sensors are deployed in order to monitor the ballistic launch site. The first three 
active ground based radar sensors (RF#1, RF#2 and RF#3 in Figure 1) are placed close 


enough to each other so as to evaluate the performance of the likelihood ratio theory 


XV 


developed in this thesis with regard to the track scoring issues. The other four active 
ground based radar sensors (RF#4 to RF#7) are placed at different random positions as 


depicted in Figure 1. 


“Geile” 


Pointer 42°54:34 26" N 13622910 91° So Eye alt 3314.03 mi 





Figure 1. Geographical area of the scenario (From:2007 Google/2007 NASA/2007 Europa 
Technologies Image) 


For each active and passive sensors located along the missile trajectory, the 
signal-to-noise ratio (SNR) is developed as a figure of merit to measure their ability to 
track the threat during its boost phase trajectory. The two IR sensors and two LIDAR 
sensors are placed as shown in Figure 1. On both satellites (SA T#¥1 and SAT#2) there is 
one passive IR sensor and one LIDAR. Figure 2 depicts the SNR for the first infrared 
sensor on SAT¥1 as a function of time. The green line represents the upper stage, which 
the missile follows after the separation of the booster. The red line represents the lower 


stage of the booster. 
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Figure 2. SNR versus Time for the First IR sensor on SAT#1. 


The use of passive infrared sensors placed on LEO satellites is very important as 
this provides the capability of early warning of the launch of a ballistic missile based on 
its radiation during the first few seconds of the flight as well as the discrimination of the 
warhead from other objects, which separate themselves from the body of the missile 
during the flight. Such objects are the booster and debris. The drawback of IR sensors is 
that they provide only angular information about the target. For this reason, an active 
laser sensor is placed on the same platform with the IR sensor in order to provide to the 


IR system the third dimension of range, enabling the ability to perform 3D data fusion. 


The calculation of the track scoring function is to identify the sensor with the best 
track file. A track score is calculated for each sensor based on the kinematics of the 
missile flight parameters and the SNR at the sensor. By using likelihood ratios, the 
optimum track file of the threat can be determined and the corresponding track file can be 
transmitted to the battle manager control in order to lead the interceptor vehicle against 
the threat. Using the optimum track file scoring techniques developed in this thesis, the 
best track file can be sent to the interceptor to guide it to the ballistic threat. This leads to 


a faster response management where the threat can be destroyed inside the territory of the 


XVil 


country that launched the missile before the missile employs any countermeasures. Figure 
3 shows the track score for three RF sensors as a function of time (or scans) and 
demonstrates that during the threat engagement, some sensors have a higher track score at 
certain times while others have a higher track scores at other positions of the threat flight. 
The objective of the track score processing is to calculate the track score for each track 


and select the optimum track that should be used to drive the interceptor. 
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Figure 3. Track score for the first three active sensors 


Based on the result of the data fusion, it is inferred that at the beginning of the 
flight the third radar (RF#3), which is closer to the threat at the beginning, has the best 
track score. As the flight time continues, the first radar (RF#1) then provides the best 
track score for the next fifty seconds. Finally, the second radar (RF#2) dominates, 


presenting the best track score until the end of the boost phase. 
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I. INTRODUCTION 


Over a decade ago, after the fall of the Soviet Union, when the world strived to 
understand why there was so much hatred among the nations, Henry Kissinger made 
these comments: 

“The end of the Cold War has made such a strategy largely irrelevant. Barely 
plausible when there was only one strategic opponent, the theory makes no sense in a 
multipolar world of proliferating nuclear powers. Mutual destruction is not likely to work 
against religious fanatics; desperate leaders may blackmail with nuclear weapons; 
blackmail or accidents could run out of control. And when these dangers materialize, the 
refusal to have made timely provisions will shake confidence in all institutions of 
government. At a minimum, the rudiments of a defense system capable of rapid 
expansion should be put into place.” [1]. 

With this, the experience of these years in our minds, the picture of the world 
right now is exactly depicted by these comments. Local conflicts among the nations have 
arisen in order to attain goals, which had been hidden for many years under the pressure 


of a bipolar control system of two super powers. 


A TRACK SCORE PROCESSING OF MULTIPLE DISSIMILAR SENSORS 


Ballistic missiles are one of the diplomatic tools used by countries to exert 
pressure on others, and fortunately, until now, their use has been restricted. What happens 
if this situation does not remain under control? At this moment, several countries around 
the world have the capability to build such missiles and exert their supremacy. Moreover, 
many of them have the capability of incorporating mass destruction warheads into these 
weapons, making the situation a grave concern, not only critical for the opponent’s rival, 
but also for all the other nations due to the possibility of a nuclear threat or accident. 

Against that threat, the most technologically-advanced, wealthy nations are 
developing advanced missile capabilities and, in parallel, are working on systems that 
provide the capability to defend their nation and their allies. The creation of an early 


warning network of different types of sensors is the goal in this regard in order to cover 


all possibilities of attack from every edge of the earth. Specifically, these networks will 
include radars, infrared (IR) sensors and light detection and ranging sensors (LIDAR) 
which will be placed on terrestrial, airborne or spaceborne platforms. 

The ballistic missile’s flight can be divided into three phases: boost, midcourse 
and terminal. The boost phase extends over the first four minutes (approximately) of the 
missile’s flight and varies according to whether or not the fuel that is used is either liquid 
or solid. The main goal of these sensors is to form a surveillance network to detect and 
track a ballistic missile threat from the launch and throughout the boost phase. The 
interception of the missile during the boost phase provides several significant advantages 
for the defender. The destruction of the missile during the boost phase provides a greater 
chance that debris and remains of the missile will fall into the territory of the country that 
launched the missile. Second, the missiles do not have the chance to deploy 
countermeasures against an interceptor. Finally, if the effort during the boost stage is not 
successful there is time to defeat the threat in the next stages. This thesis work is focused 


entirely on the boost phase of the threat. 


B. PREVIOUS WORK AND OBJECTIVE 


This thesis is a continuation of the work started by Rakdham [2] who examined a 
scenario involving the launch of ballistic missiles against the United States from North 
Korea, developed the multiple hypothesis tracking and extended the theory of the linear 
assignment problem (LAP). The LAP theory was initially developed by Danchick and 
Newman [3] to provide an alternate approach to solving the related data association 
problem by minimizing computational requirements. Rakdham [2] also used the 
IMPLUSE model from NASIC in his work. 

The main goal of this thesis is to examine the issue of track scoring for the 
different dissimilar sensors within the surveillance network and the data fusion process of 
the raw data being delivered over the surveillance network. Data fusion is then the 
selection and the compromise of the accumulated data, which will be transferred to the 
battle manager control (BMC), turning the raw data being collected into decision level 


knowledge for the interceptor guidance. 


C. PRINCIPAL CONTRIBUTION 


The major tools employed in this effort include the IMPULSE simulation tool for 
launch of the ballistic threats and the MATLAB simulation software. The IMPULSE 
simulation tool was developed by the National Air Space Intelligence Center (NASIC) 
and provides the trajectory for a given ballistic missile, which is necessary for the 
creation of the required data files. These data files are then processed by the MATLAB 
code developed in [2], which (a) modeled each of the RF sensors, (b) performed the 
multiple hypothesis tracking and (c) developed the linear assignment problem and 
solution. This thesis extends the effort in [2] by including additional IR and LIDAR 
sensor systems within the surveillance architecture and also increasing the number of RF 


sensors that are used. 


The IMPULSE program was used as a tool for analyzing the trajectory file of 
ballistic missiles. This program has the capability to provide all aspects of the kinematics 
during the flight of the missile, such as the position, the velocity, the acceleration, the 
change of the mass and the change of the drag force on the missile. Also it provides the 


capability of virtually viewing the flight as a function of time. 


Next, the radar cross section (RCS) of the ballistic missile was examined in order 
to deploy a more dynamic aspect of the radar equation, which would be necessary in the 
tracking function and data fusion process. The calculation of the RCS is based on the 


geometry of the trajectory and the incident angle of the beam on the missile. 


After this step, more sensors are added in the vicinity of the launch site in order to 
evaluate their response against the ballistic threat missile. Two infrared (IR) and two 
LIDAR sensors are placed on two spaceborne platforms (low orbit satellites). That is, 
every spaceborne platform contains one IR and one LIDAR sensor. The satellites are 
deployed close to the vicinity of the launched missile to provide good coverage. Finally, 
having the responses of each sensor within the surveillance network, the track scoring 


and data fusion signal processing could be developed. 


D. THESIS OUTLINE 


Chapter II introduces the sensors for the data fusion system developed in this 
thesis and also discusses their implementation in the MATLAB code. It begins with a 
definition of the SNR of the IR sensors and analyzes all the parameters that must be taken 
into account. It also develops the IR sensor model, which will be placed on the low-earth 
orbit (LEO) satellite. The second section of this chapter analyzes the LIDAR sensor SNR. 
At the beginning, the SNR formula of the LIDAR sensors is developed by covering all 
the necessary aspects. After the definition of the SNR, the LIDAR model of the 
simulation is introduced. The chapter ends by providing the parameters of the deployed 


radars and the evaluation of the missile’s RCS. 


Chapter III discusses the data fusion of the tracks of deployed sensors. The theory 
of likelihood ratio is developed based on the contribution of the kinematics and the sensor 


signals. 


Chapter IV depicts the scenario considered. Specifically it provides the positions 
of the sensors and their responses to the ballistic missile. Also, the result of the data 


fusion problem for the radar sensors is presented. 
Chapter V presents a summary and the concluding remarks. 


Finally, Appendix presents the flow chart of the MATLAB code used in the 


simulation. 


Il. SENSOR DESIGN 


This chapter presents two space based sensors (infrared and light detection and 
ranging system) and one radio frenquency sensor (radar). The signal to noise ratio and 


design considerations of these sensors are discussed in depth. 
A. IR SENSOR 


The signal-to-noise ratio (SNR) is the figure of merit used in this work in order to 
evaluate the track quality of an infrared sensor. Also, it is known that the SNR is the 
basic figure of merit for all the sensors with respect to the design and evaluation of their 
ability to perform target tracking and detection. In order to define the SNR of a system, it 


is necessary to determine the ratio of the signal power to the noise power. 
1. SNR of an IR Sensor 


First, SNR is defined for the best case in which the background is not taken under 
consideration, assuming a blue sky background. The signal power which a sensor can 


calculate from a point source is given by: 





H= = Wem” (2.1) 


where J is the radiation intensity in Wsr?! and R is the range to the target in centimeters. 
y g g 


Note that (2.1) does not consider any atmospheric absorption. Including the 
transmission losses due to atmospheric absorptions in the (2.1) yields: 
J 
R 





H =7(R) W cm” (2.2) 
The effect of transmission loss in an IR system is a very important consideration 
in determining the performance of such a system. There are three major parameters that 


characterize atmospheric attenuation [4]: space extinction, emission, and turbulence. 


The transmittance is a measure of the atmospheric attenuation of the energy from 


a light source or target and is given by [3] 


(R)=1-€ (2.3) 


where R is the range, and / is the sum of the two attenuation parameters: absorption 7 


a? 


and scattering #.. The transmittance though the atmosphere is presented in Figure 4. 
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Figure 4. Atmospheric absorption and scattering (From: 10) 


As shown in Figure 4, there are three IR bands available for applications: short- 


wave IR (SWIR), which is in the 1-2.7 44m wavelength, the midwave IR (MWIR), which 
is in the 3-5 42m wavelength, and the longwave IR (LWIR), which is in the 8-14 wm 
wavelength. For each one of these bands, the parameter ( is very complicated to 
determine with respect to all the atmospheric characteristics, such as rain, fog, clouds, 


haze, and so on. Standard computer models for obtaining parameter @ are available 


[4],[5]. 


The amount of energy that an IR sensor can receive from a source depends on the 
aperture of its lens. Given that the aperture is d (m), the amount of energy that the sensor 
receives is called the irradiance and is given by: 

J nd’ 


— [W] (2.4) 


S =7(R)— 
(Re 4 





As with all sensor systems, such as radar, the energy measured by the system at 
the antenna is not the final energy the system uses to calculate target position. This is 
because there are many sources of noise inside the system. In fact, in an IR system, there 
are two types of noise: the background noise, which it is not considered initially for the 


calculation of the SNR, and the noise in the sensor. 


The noise in the sensor can be categorized as Johnson noise, shot noise, 
generation-recombination noise, temperature noise, microphonic noise, 1/f noise, and 
finally electronic-interface noise [6]. 

Then the internally generated noise power, called noise equivalent power (NEP), 
is given by [6]: 

NEP =1/D (2.5) 


where D_ is the Detectivity of the sensor. Note that this detectivity can also be a 


function of wavelength D(A) . 


Finally, the specific Detectivity, D°, another figure of merit for IR systems is 


defined as [3]: 


_(AAfY” 
NEP 


p* (2.6) 


where Af is the equivalent noise bandwidth, and A, is the effective sensor area. 


All of the system losses, such as losses due to optics, electronics and detector, are 
accumulated in a loss factor L, which is applied to (2.4) and results in the signal power at 


the detector as 


2 
S =r(R) i md 
R 4 








L [W] (2.7) 


In the following analysis, the sensor is considered to be operating in a staring 
model. The advantage of this model is that all of the detectors are considered to cover the 
field of view (FOV) for the entire time frame [6]. In this case, the number of pixels is 


equal to the number of detectors in the array, so the time frame (7,) to scan the entire 


7 


image is equal to the dwell time (z | , the time required to read all of the pixels, which is 
expressed by 
T, =T,. 
(2.8) 
As a result, the equivalent noise bandwidth Af , which is also the signal bandwidth 
[6] page 489, is equal to: 


1 1 


The instant field of view (IFOV) @, can be obtained from the detector area and 


the focal number( f /#)d and is given by: 


ae: ieee (2.10) 


[ray 
where A, is the detector area. Finally, combining (2.6) to (2.10) allows us to express the 


SNR as 





(2.11) 


S — ndID*L (nt, ) 
N 4(fA)R?\ Q 


where Q is the solid angle, which the FPA observes, in steradians. 
Zi Space-borne IR Sensor Design 


The simulation developed in this thesis employs two IR spaceborne sensors, 
located at different positions around the Earth, both in the vicinity of the launched 
missile. The satellites are in a low-earth, sun-synchronous orbit (LEO) at an altitude of 
705 km, crossing the equator at the same local time on each pass. Their ground velocity is 
6.76 km/sec; therefore, it takes approximately 90 min for a single revolution around the 
Earth and requires 16 days to revisit the same location on the Earth’s surface. Note that 


the satellites are assumed to be stationary for this simulation. Otherwise, the simulation 
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has to take into account the relative motion between the spaceborne platforms and the 
ballistic missile, which makes the problem more complicated and is beyond the scope of 
this work. The modeling parameters used for the two IR spaceborne sensor satellites are 
taken from the Landsat 7 satellite [7]. This is National Aeronautics and Space 
Administration’s (NASA) most recent model and is used for commercial reasons, 
assuming that a rotating telescope with 1000x1000 detectors in the focal plane array for 
better resolution of the target, which means better discrimination of multi targets at the 
same area. Another assumption is that only band number 6 (out of eight bands) in SWIR 


is used. The Landsat 7 parameters are summarized in the Table 1 [7]. 









































Parameters Units 
Ground resolution 185m 

Wavelength 2-2.35 um (SWIR) 
Spectral bandwidth 10 nm 

IFOV 262.4 pt radians 

Telescope diameter 55mm (2in) 

Optical f/# 6 

Number of detector per spectral band 1000 x 1000 
Detectivity 3.5x 10'° cm-Hz'?/W 
Loss factor 0.085 
Table 1. | Sensors Parameters 


The basic source of the infrared energy from a ballistic missile is the hot exhaust 
plume. The other sources of infrared energy are the exhaust nozzle of the engine and the 
aerodynamic heating of the missile [8]. This work does not consider the latter two 


parameters, only focusing on the radiation from missile’s plume. An estimate of the 


magnitude of this radiation intensity can be made by taking into account the size and the 


temperature of the associated plume observed by the detector [8]. 


Right after the launch, the size of the plume is 4 m in diameter and about 50 m 
long, and the temperature of the plume at the nozzle is about 1800 K. Using an average 
temperature of 1400 K, the highest radiance exitance can be calculated from the 
blackbody equation and is 7x 104 Wm?-pm’' at a wavelength A of 2.1 um. Assuming that 
the plume is a blackbody and radiates isotropically, then the resultant radiation intensity 
is approximately 3 MW/sr. As the target rises, it comes closer to the sensor and results in 
an increase in the radiation intensity due to an expansion in the size of what the detector 
can observe as well as the associated decrease in the atmospheric absorption. For the 
purposes of this simulation, the radiation intensity is assumed to be a constant at its peak 


value of 3 MW/sr. 


As mentioned in the beginning of the section, the detector is a rotating telescope 
that uses the focal plain array (FPA) technology. The detector is a cooled mercury 
cadmium telluride (HgCdTe). FPA technology is very important for the discrimination of 
ballistic missile at the boost face, using more color bands [9]. Specifically, two bands are 
used, one before and one after the burnout. The HgCdTe provides a detection capability 


across many wavelengths from the SWIR to very long wavelength (14 um). 


The detectivity of the IR sensors used in this simulation is estimated from Figure 
5 for a wavelength of 2.1 um. The values of detectivity in this figure were derived from 
an analysis of commercial detectors, under background- limited conditions, at the 
temperature of 300 K and an IFOV of 180° [10]. In this case, it is also assumed that a 
filtered, cooled detector at a temperature of T = 77 K was used. Taking into account all 


these assumptions, the detectivity from Figure 5 is equal to 1.5 10'° em-Hz'”/W. 
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Figure 5. D* for commercial detectors (From: [10]) 


At this point, the energy from the terrain, called background IR, is analyzed. The 
analysis covers only the first four minutes of the missiles’ trajectories in the boost-face. 
During this interval, the background energy comes only from the terrain as the missile 
has not reached the peak of its altitude; therefore, the satellite looks only at the terrain. 
Otherwise, as the missile rises into the space, the background will change making the 


problem more complicated. 


At the wavelengths of 2 to 2.35 um, the background energy from the terrain is 
generated by the reflection of the solar radiation [11]. This energy is calculated assuming 
that the terrain is a blackbody source at a temperature of T = 300 K. For wavelengths 
between 3-5 um, the background energy is considered negligible due to the atmospheric 
absorption [11]. It takes large values of background energy for wavelength above the 5 


um, which does not affect this work. Figure 6 depicts the energy from the background. 
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Figure 6. IR signal from the missile and the background (After: [11]) 


The background energy from the terrain at a wavelength of 2.1 um is calculated 


as: 


he? 1 





=5220 Wem? -um (2.12) 


he 


ekT —] 


As a matter of fact, this is not the background energy, which the detector 


measures because the detector is cooled. Figure 7 depicts the geometry of a cold detector. 


12 









Detector background 





Background 
Extended 
source 











Cold aperture 


Figure 7. Cold detector geometry ( From:[8]) 


The detector background is then given by 
= 7M (A)Q = 30.33 dB (2.13) 


E, background 
The total solid angle coverage is computed by the elevation coverage angle and 


the azimuth coverage. The elevation coverage angle: 


€= arctan = =14.7° (2.14) 
705 


The azimuth coverage angle is given by: 


h= arctan ane =14.7° (2.15) 
705 


Finally, the total solid angle : 


S29 0658 a (2.16) 


732 


O= 





The scan geometry of the employed spaceborne sensors is shown in Figure 8. 
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Figure 8. Scan geometry (After: [10]) 
The number of the resolution cells (pixels) is calculated as: 


Q 0.0658 


P=—~ =—————_.,, ~ 1 Mpixels (2.17) 


ay (262.4-10°) 


Figure 9 depicts the SNR value of the IR sensor #1 placed at the satellite platform 


SAT#1 by using (2.11) for SNR in MATLAB. In order to calculate the SNR, it is 
assumed that the range of the missile is provided to the IR sensor by the laser sensor as 


described in the following section. The simulation calculates the SNR for a 3D geometry 


instead of a 2D geometry. 
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Figure 9. SNR versus time for a IR system at the platform SAT#1 


From Figure 9, the SNR is about 25 dB at the beginning of the flight and is a 
reasonable value for detecting a ballistic missile. Also, it can be mentioned that after the 
separation from the booster at time f=65sec, the missile rises closer to the sensor 


having a peak SNR value of 38 dB. 
B. LIGHT DETECTION AND RANGING SYSTEM 


This section presents a derivation of the SNR of the light detection and ranging 


systems (LIDAR). Also, spaceborne LIDAR design issues are briefly discussed. 
1. SNR of a LIDAR Sensor 


LIDAR systems are similar to radar systems. They are active systems as they 
illuminate the target in order to achieve the parameters of the target position. Compared 


to a conventional radar, the beamwidth of a LIDAR is very small due to its operation at 
i 


very high frequencies. This trait gives them the capability of high resolution. On the other 
hand, due to the narrow beamwidth, they have limited search capability. For this reason, 
they are used in conjunction with other electro-optical systems. An example of these 
systems is the laser range finder (LRF) which is used together with a forward looking 
infra-red (FLIR). LRF provides the range, in addition to the angular data from the FLIR, 


to a Kalman filter to form a 3D track. 


In the design of a LIDAR system, the atmospheric absorption has to be taken into 
account. As shown in Figure 2, the electromagnetic radiation through the atmosphere at 
um wavelengths specifies two bands where the laser radar can be applied. These 
windows are at the 1.06 um and 10.6 um wavelengths. At 1.06 tm, the system uses 
neodymium YAG crystal lasers, and at the 10.6 um the system uses CO) gas lasers. 


The following derivation for the SNR is taken from [10]. To define the SNR, it is 
necessary to calculate the signal power received by a LIDAR detector. The distribution of 
signal power over a spherical area is given by: 


P 
An R? 





(2.18) 


where P. is the transmitted power and R is the distance of the target. 


t 


The focusing effect of the lens in order to direct the transmitted power in a 
specified direction is given by: 
Ar 


T 72 
— Ony 


4 





L= (2.19) 


where @;,, is the beamwidth. 


In this case, it is assumed that the ballistic missile is a Lambertian target, so the 
power reflected by the target is distributed over a hemispherical area. For this reason, the 


signal power is divided by the factor 1/nR’. 
Finally the signal power is given by: 


PAT, 
S= 5 
27aR 





(2.20) 


16 


where p is the target backscattering coefficient, 7, is the optical efficiency and A is the 
lens aperture area. 
The above form covers the case in which the laser beam is much smaller than the 


target area. Otherwise, we should take the target cross section o into account. Then the 


signal power takes the form: 


2PAot 
Se (2.21) 
mR Onw 
where 0,,,1S the parameters of lens. By substituting@,,,=2/D, where/is the 


wavelength and D is the lens diameter, (2.21) takes the form: 


> 8PA°oT,0 


: wR)? 


(2.22) 


From [10], the detector noise level is determined. The LIDAR detectors are of 
two types. The direct (or nonconherent) sensors, which accumulate the energy in a way 
similar to the electro-optic sensor, and the coherent (heterodyne) sensor in which a phase 
difference between the transmitted and the received signals is utilized. In this case, the 
range is measured by the Doppler phenomenon. In this simulation, coherent sensors are 


used. 


The NEP for a LIDAR is given by the form: 


NEP =(A,Af)'? /D* (2.23) 


where D* is the detectivity, Af is the receiver bandwidth , and A, is the detector area . 


Finally, by taking into account (2.16) and the noise level for a heterodyne receiver 


as given by [10], NEP is obtained as: 


NEP = (hf / j)\(B) (2.24) 
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where j is the quantum efficiency, B is the receiver bandwidth, the Planck’s constant 


and f is frequency . 


The SNR of the LIDAR system is then given by: 


s 8PApt,o 


SNR = = 
NEP w°R*(hf / j)(B) 





(2.25) 


As mentioned at the beginning of Chapter II, this value will be the figure of merit 


in order to evaluate each track of the spaceborne sensors. 
2. Spaceborne LIDAR Sensor Design 


This section discusses the design of spaceborne laser sensors. An NdYAD laser 
system, operating at a wavelength of 1.06um, is placed on a LEO satellite at an altitude 
of 705 km around the Earth (along with an IR sensor). The parameters of a typical sensor 


of this type are listed in Table 2 [10, chapter 8]. 





























Parameters Units 
Detector area 6.25 x 10% cm? 
Optics diameter 20 cm 

Wavelength 1.06 um 

Pulse power 2 MW 

Pulsewidth 10 ns 

Target reflectivity 0.1 
Detectivity D* 10'° cm-Hz'?/W 
f lit 2 
To: optical efficiency 0.5 














Table 2. | Parameters of LIDAR system[Ref 10, ch 8] 
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The inverse of the pulse width defines the receiver bandwidth B =1/r =10° Hz. 

From (2.22), noise equivalent power of the receiver can be written as: 
NEP = [(6.25x 10“)(10°)]!7/10''=2.5 x 10°W (2.26) 
The range of the missile is calculated by the MATLAB function, 
LIDARobservation, by using the inputs from the IMPULSE simulation program. In the 


simulation, it is assumed that the ballistic missile with dimensions of 20 m length and 1.2 


m of diameter is a perfect cylinder. As a result, RCS iso = 32 dBsm. 


In Figure 10, the SNR is depicted as a function of time for the LIDAR sensor #1 
placed on platform SAT#1 of the simulation. 
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Figure 10. SNR versus time for the LIDAR sensor on SAT #1 


C. RADAR DESIGN 


Seven radars, which are placed at different locations relative to the launch 
position, use the same configuration in order to make the analysis easier to verify. In this 


section the radar model originally developed in [2] and used in the simulation is reviewed 
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with a summary of the parameters given below in Table 3. The single pulse, signal to 
noise ratio (SNR), for each radar system can be expressed as: 


PG?’ 
Gi ees 


2 227 
(47) kT, BR* 220) 


where 7.=T,+T,, T, is the receiver noise temperature and 7, is the antenna noise 


temperature which in the work reported here assumed to be 290 K as was the assumption 
in [2].P is the peak transmitted power, G, is the antenna gain, B is the bandwidth of 
the receiver,o is the radar cross section of the ballistic threat target, k is Boltzmann’s 
constant. Due to the large unambiguous range required, the pulse repetition frequency 


(PRF) is modeled as PRF=150Hz. 


Table 3 lists the parameters of the radar considered in the simulation. 














Parameters Values 
Transmitted Power 1 MW 
Carrier Frequency 10 GHz 
Antenna Gain 42 dB 





Radar operating Bandwidth | 15 MHZ 
Receiver Noise Temperature | 290 K 
Table 3. Radars parameters [2] 

















The radar cross section (RCS) of the missile during the flight is calculated by the 
program POFacets [12]. The missile used in the simulation is a sample unclassified 
model. To determine RCS of the missile, it is assumed that the missile is a circular 


cylinder with a radius of 0.60 cm and a length of 20 m. 


The angle 6 between missile’s trajectory vector and the distance vector from the 
radar to missile, which is required in order to define the RCS function, is determined by 
the MATLAB function RFobseravation by using the geometry of the Figure 11. From 
Figure 11, it is inferred that the angle of incidence is computed as the inner product 
between the vector of distance from the sensor to missile and the vector of the curve. The 
vector of the curve is defined from the difference of the successive points of the 
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trajectory. The RCS determination is a 3D problem. In this simulation, it is assumed that 
the change of the angle takes place only on the x-y axes without taking the third 


dimension into account. 


weet tees 
PT Loh alata telat dS) 
Phd . 
Po 
oe? 
. 
* 
.- 
. 


ga ieeeimctaren <<*** curve of the trajectory 
Target (ballistic missile) ° 


Rs LOS 
.” Incident angle 


R: distance between target and radar 


‘2 Radar 


Figure 11. Calculation of incident angle 0 


Using the above geometry, the MATLAB function RFobservation calculates the 


incident angle 6 for the first radar sensor and is shown in Figure 12. 
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Figure 12. Calculation of incident angle 0 


Having determined the incident angle and using the POFacets [12] code, the RCS 
of the ballistic missile is obtained as shown in Figures 13 and 14. 

Figure 13 is a polar plot of the RCS as a function of the angles of incidence. In the 
vicinity of angle @=90°, which has been calculated from the previous analysis, the 


values of the RCS are around 36 dBsm. 
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Figure 13. Polar plot of the Ballistic Missile RCS 
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Figure 14. Linear plot of the Ballistic Missile RCS 


Figure 14 is a plot of the ballistic missile RCS with respect to the angle of 
incidence. From the above two plots, it is inferred that the RCS for an angle of incidence 
in the vicinity of 90° is equal to 36 dBsm. 
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In summary, this chapter described the parameters of the three different types of 
sensors used in the simulation. The parameters for the IR and LIDAR sensors are listed. 
Also, the parameters for seven radars at sea level are listed. Next, Chapter III describes 
the data fusion process and the theory of the likelihood ratios, which is used in the 


simulation for evaluating the performance of the radars. 
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HiIl. DATA FUSION PROCESS 


Data fusion process is one of the most critical parts of an integrated network of 
different types of sensors placed on terrestrial, airborne, or spaceborne platforms. The 
other processes that are part of the data fusion process are the association of the targets, 
feature extraction and identity declaration. Specifically, data fusion is the procedure of 
collection, combination and evaluation of information from different sources (sensors), 
seeking to estimate and predict the level of threat against our friendly positions and units. 
There are many techniques of collecting and compromising data which depends on the 
purpose and the design of the system. All these techniques are basically similar as they 
follow the same philosophy. The difference each time is the place where the association 
and data fusion take place in an integrated system. In this simulation the series of the 


procedure is described in Figure 15. [13] 
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Figure 15. Procedure of identity declaration (From: 13) 


From Figure 15, each sensor makes its own association concerning the multiple 
targets that are present. The association processing of the targets used in this work is 
according to [2]. Continuing on, the results of the association of each sensor are sent to a 


fusion center where the data fusion takes place. Then, based on the result of the data 
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fusion, the feature extraction and finally the declaration of the identity of the target is 
derived. With other techniques the data fusion comes first in the center and is responsible 
for the association of the targets as well as the feature extraction and the identity 


declaration. 


There are many advantages to fusion of data from multiple, dissimilar sensors 
placed in different locations. First of all, the combination of different observations results 
in an improved estimation of the target’s trajectory (velocity, position, acceleration). That 
means, statistically, using the results of N sensors, the estimated position of the target 
can be improved by a factor N'” [14]. Also, by using observations from many sensors, an 
improved estimated position of the target can be achieved. As an example, consider the 
triangulation between two sensors that observe the same target [13]. Second, by using 
different types of sensors, the fusion performance is improved as the information of one 
sensor can be used to overcome the deficiencies of other types of sensors. For example, 
the case of an IR sensor alerting a radar to confirm the presence of a target can improve 
the probability of detection [4, Chapter15]. The IR provides the detection of the target 
and the angular direction of the target to the radar. Then the radar, by using a longer 
standard coherent waveform, seeks to confirm the target and attain the range as well as 
the range rate data. The radar, based on the SNR of the IR system, can lower its threshold 
in order to adjust the detection performance on the target. As a result, not only an 
improved detection of the targets takes place, but also the attributes of each sensor are 
utilized. This is an example of the interpolation techniques, which can take place between 
different types of sensor (FLIR, LIDAR, airborne sensor, space board systems), in order 


to cover any aspect of the problem. 


Figure 16 depicts the sequence of the procedure of the interpolation techniques 
described above using a radar and an infrared search and track (IRST) system. In this 
example, the radar and the IRST use inter-sensor cueing and inter-sensor data sharing to 
confirm the presence of the target. The data fusion process develops a central tracking 
file, which is supplied to the sensor management processor. The sensor management 
process then feeds the processed target results back to the radar and the IRST for 
refinement. 


26 


RADAR 





Signal 
r processing 











Tracking 


Observation 





Confirmed 


















































sensor data 
Central 

: level 
I Tracking 
, data 
a a all geeiet eal Data association | Data fusion ——> 
; JN WV 
‘ ; Confirmed , 
I oben eee Tracking sensor data . 
j Processing Processing 
I ¥ 
| IRST Sensor 
iat ta faa ho i eo re tt ee ee ereceniens 

Figure 16. The role of coupled fusion (After: [4]) 

A. TRACK SCORE FUNCTION 


In this section, the likelihood ratio of the kinematic parameters, and the signal 


parameters, which contribute during the measurement process, are developed. The 


likelihood ratios are the figure of merit for this work in order to evaluate the track 


reliability. For this reason, it is called the track score function. 


The likelihood ratio for a combination of data, taking the priori probability data 


into account, is given by [15, page 73]: 


A(H 


H,) 


= p(D/H,)Po(H,) A FE, 





1? 


ai 


p(D/H,)Po(H,)  P. 


(3.1) 


where H,is the true target hypothesis, H,is the false target hypothesis, P. is the 
probability of true target, P, is the probability of false target, p(D/#H,) is the probability 


density function assuming that the H, is correct, Po(H,) is priori probability of H,. 


The performance of the tracking system should be independent of the behavior of 
the target. The system is evaluated for its ability to respond instantly to any change of the 


trajectory of the target either in velocity or the direction. 


Taking K scans of data, where the measurement error for one is not related to the 
error from the previous scan, the likelihood ratio (A) can be given as the product of the 
likelihood ratio due to the kinematics (A,) and the likelihood ratio due to the signal (As) 


is given by [4, Chapter 15]: 


K 
A(K)=Ao IT A,A, (3.2) 
where 
Po(H,) 


The initial track score is dependent on the first result of the measurement of the 
system. As the system starts to count from zero, the kinematics parameters do not play 
any role in the initial track score. The initial track score can be expressed as [4, Chapter 
6]: 


_ Po(H,)P(Det/H,) p(y, / Det, H,) 


ADS Po(H,)P(Det/ H,) p(y, / Det, H,) 





(3.3) 


where y, is the total signal related data (in this work, the signal-related data will be the 
SNR of the target), p(y,/Det,H,)is the probability density function of detection 
assuming that H, is correct, p(y,/Det,H,)is the probability density function of 
detection assuming that H,, is false alarm, and Po(H,) is the initial probability of target’s 


presence, which, is given by: 
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Po(H,) = ByrVe (3.4) 
where V, is the measurement volume (we assume that every true target detection or 
false alarm takes part inside of this volume; in this work, V. is equal to one), Bex is a 


new target density (assumed equal to unity). 


The case of false alarm is given by the P,(H,) =1. Additionally, P,= P(Det/H,) 


is the probability of detection and P,,= P(Det/ H,)is the probability of false alarm. 


Combining (3.3) and (3.4) and taking P,(H,) into account, the initial track score 


is given by: 





A(l) = ByrVePy p(y, / Det, H,) 








(3.5) 
Pa PCy, / Det, Hy) 
Taking the logarithm of (3.5) yields: 
In A(1) = Inf By, VJ + Inf [np 2O Pea) (3.6) 
Pra p(y, / Det, Hy) 


In this work, the initial track score is assumed to be A(1)=0 . 


1. Kinematics Related Likelihood Ratio 


To determine track score based on the kinematic parameters, it is necessary to 
introduce gating. Gating is a method of determining which observations are the best to 
update the tracks. Gating defines cells where the measurements and the observations are 
included in order to find the best match among them. The scope of gating is to define as 
small an area as possible, which will include both the measurements and the observation, 
minimizing the possibility of error. There are many techniques of gating. The cell could 


have a rectanguloid or an ellipsoid. 
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The term global nearest neighbor (GNN) is the method used in this work. GNN is 
the method of association in which each measurement in a track is associated with the 
nearest neighbor in order to achieve the best update by its closest observation [4, Chapter 


1]. Figure 17 depicts the problem of gating. 


Gate Gate 





1,2,3: Observation positions 


P,P,: Predicted target position 


ae: Distance from P, to 2 


Figure 17. Gating (From[4}]) 


The distance d*y) from predicted target position P,to 2, called normalized 


statistical distance, is given by [4, Chapter 6]: 


dy =y' S'y (3.7) 
where yis innovation (the difference between the actual and expected measurements at 


every track), S is measurement residual covariance matrix. The innovation y is given by 


30 


y= y(k)—AC(k/ kD) (3.8) 
where / is the nonlinear measurement function for the case of the extended Kalman 
filter, x is the state vector, X is the state estimate vector, y(k) = h(x(k))+v(k), and v is 
white Gaussian noise with covariance matrix R and mean zero.. 

In a nonlinear function like the extended Kalman filter, S~' is given by the 
following form [4]: 

S"'=HPH'+R (3.9) 
where P is the prediction covariance matrix and H is the measurement matrix. 

Taking the gating into account and assuming a Gaussian distribution for the target 
returns as well as a uniform distribution for the false alarms, the likelihood ratio (A,) for 
the kinematic part is defined as: 

_ V, el !2) 


Wi pee oe 
k (nyt 


(3.10) 


where m is the measurement dimension (in this work, it is set to one) and |S | is the 


determinant of the measurement residual covariance matrix S$ . 
2. Signal Related Likelihood Ratio 
Based on likelihood ratio (As) of the signal-related data, the track score is given 


by [4, Chapter 6]: 


_ P(Det/H,) p(y, / Det, A,) 
*  P(Det/H,) p(y, / Det, H,) 





(3.11) 


Taking both Pp, and P,, into account, (3.11) can be written as: 


rem P, p(y, / Det, H,) (3.12) 
Pp. / Det, Hy) 





S| 


In the case of the radar, when the SNR is available, the signal likelihood ratio is 
given by (3.11). The Swerling model 3 is typical for the missiles under consideration [4, 


Chapter 2]. The probability density function under H, is given by 





PS) ees -y ] (3.13) 


P(y/ A) = = 
(1+ ©/2)2 14+0/2 


where y is SNR data and © is the average of the SNR. 


Probability density function under H, is given by: 


P(y/ Hy) =exp(-y) (3.14) 


For an IR sensor, the derivation of the probability density function is more 
complicated. A constant false alarm rate (CFAR) detector is used in order to detect 
unresolved targets. CFAR detectors are designed to maintain a constant false alarm rate 
for a specific period of time. For this reason, a CFAR detector normalizes the background 
clutter by dividing the filter output by sample standard deviation, which is assumed to be 
known or can be obtained through a Monte Carlo simulation. How this occurs in practice 
is explained in [15]. Finally, it is known that, by using a sample standard deviation, the 
SNR follows a t-distribution instead of a Gaussian distribution. That means, the SNR has 
larger tails and larger values above the threshold. As a consequence, the probability 


density function under #7, can be written as 


exp(—s2(1—a?)/2) : 
I(s, 3.15 
[Nar /DOs yy Ne oor 


as y’ /N 
l+y?/N 


I(s,a) = | (z+as)%e* dz 


—as 


P(y/H,) = 





where 


32 


and NV is the number of degrees of freedom. For IR systems, N is relatively small. In 
order to calculate the maximum likelihood value of N , the number of data samples in the 
tail of the distribution should be calculated first [15]. In this work it is calculated by 


testing the system performance. 


The probability density function under H, is given by: 


1 T(N+1/2 1 ’ 
POH) = ae oa (N+1)/2 


3.16 
JNx Y(N/2) oy 


In summary, this chapter described the theory of the likelihood ratios used in the 
simulation. The likelihood ratio formula for the kinematic and the signal-related 
parameters of the missile is defined in order to calculate the track scoring function for a 


radar. Next, Chapter IV demonstrates the simulation and the results. 
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IV. SIMULATION RESULTS 


In this chapter, the simulation of a missile tracking is presented; the IMPULSE 
model and MATLAB are the tools used to implement the simulation. Specifically, the 
results from the IMPULSE program are used by the MATLAB code from [2] in order to 
simulate the tracking problem. The simulation results of the sensor performance based on 
the likelihood estimation are presented. In this simulation, an unclassified model of the 
ballistic missile is used. A detailed description of the requirements (inputs) and 


capabilities (outputs) of the IMPPULSE program are outlined in [2]. 


Figure 18 depicts the geographical area as well as the positions of the sensors. 
The geographical area of the scenario chosen is arbitrary. The ballistic missile covers the 


distance from Asia to North America in approximately at 40 minutes. 


“"Gog@le” 





Eye alt 3314.03 mi 


Figure 18. Geographical area of the scenario (From: 2007 Google/2007 NASA/2007 
Europa Technologies Image) 
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In Figure 18, the seven active ground based radars (RF#) and satellites (SAT?) are 
placed close to the missile launch position. The seven radars have the same parameters 
(transmitted power, frequency, beamwidth, and PRF) and are at different positions in 
order to examine the tracking and fusion performance. More radars can be used along the 
trajectory to cover the entire flight of the missile. The sensor radars used in the simulation 


are adequate for covering the area depicted in Figure 18. 


The IR and LIDAR sensors are placed on the satellite platforms. It is assumed that 
there are no interactions between them, and the transfer of the data between the two 


sensors takes place instantly. 


Table 4 lists the positions of the sensors in the vicinity of the launched ballistic 
missile. The radars are placed at the sea level, but the IR and LIDAR are placed at an 
altitude of 705 km above the sea level. 
























































a/a_| Type of sensors Position 
Latitude Longitude Amplitude 
I RF #1 131°4644 E 40°50 45 N Sea level 
2 RF#H 132°46 44 E 41°40 45 N Sea level 
3 RF #3 131°00 00 E 40°00 00 N Sea level 
4 RF #4 130°00 00 E 35°0000 N Sea level 
5 RF #5 135°00 00 E 40° 00:00 N Sea level 
6 RF #6 140°00 00 E 45°(00 00 N Sea level 
7 RF #7 145°00 00 E 50°00 00 N Sea level 
8 IR/LIDAR #1 | 130°0000 E 35°3000 N 705 km 
9 IR/LIDAR #2. | 125°0000 E 45° 3000 N 705 km 
Table 4. _ Sensor positions 


Figure 19 depicts the same geography of the scenario as in Figure 18 but from a 
different angle of view. In this depiction, the altitude of the satellites and their place 


relative to the launch position, L, is clearly presented. 
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Figure 19. Geographical area of the scenario from a different angle of view (From: 
2007 Google/2007 NASA/2007 Europa Technologies Image) 


Figure 20 depicts the trajectory of the missile as generated by the IMPULSE 
program. The type of the missile is a “Sample Unclassified Model.” The flight of the 
missile takes about 26 minutes to cover the distance from Asia to the North American 
continent. Of the two red lines in Figure 20, one represents the actual trajectory in the air 
and the second is the footprint on the earth. This graph covers the first 402 sec of the 
flight. The separation of the booster takes place approximately 65 sec after the launch 
(when the blue line becomes red) and its trajectory after the separation is the red line. The 


missile launch position is N40° 51’ 17’’ and E 129 ® 39’58”’ moving toward the East [2]. 
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av ‘gate Mode Simulation time: 402.080 Sec. 
Selected: Selected: None 


_ Actuals 
Sac jectory 


Footprint 





Figure 20. Trajectory of the Ballistic missile 


The parameters of the launched missile selected in the IMPULSE program are 
given in Table 5 [2]. 


























Parameters Inputs 
Type ‘SampleUnclassModels’ 
Missile “(U) BOOST Unclassified Sample” 
Latitude N40° 51? 17” 
Longitude E 129° 39°58”’ 
Altitude 20m 
Azimuth 23° 
Kick angle Li? 














Table 5. — Ballistic Missile parameters. IMPULSE program inputs. 
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Figure 21 depicts a panoramic view of the earth and the trajectory of the missile 


as well as its footprint on the ground for an interval of 1239 sec. 


Navigate Mode ME Simulation time: 1239.00 Sec. 


Selected: ° eae 





Figure 21. Panoramic view of the ballistic missile’s trajectory. 


A. RADAR PERFORMANCE 


In this section, the results of individual radars as well as the results of data fusion 


based on the track score mechanism developed in Chapter III are reported. 


Figure 22 shows the trajectory of the missile as tracked by RF #1. The red spot on 
this 3D graph is the position of the radar relative to the missile launch site. The red line 
on the trajectory is the output of the IMPULSE program. The blue line represents the 
tracking of the radar. The radar follows the trajectory of the missile closely. At 65 sec, 
the separation of the booster takes place, which is depicted on the graph by the small 


curved line. 
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Figure 22. 


Vertical 


The trajectories from the other six radars deployed in this simulation are depicted 


in Figures 23 and 24. 
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Figure 23. Ballistic missile trajectories from Radars #2, #3, and #4 
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Figure 24. Ballistic missile trajectories from radars #5, #6, and #7. 
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After the separation of the booster at t =65 sec, the radar follows the missile and 
the booster. The missile’s trajectory (green dashed line) is the upper stage. Booster’s 
trajectory is (red dashed line) is the lower stage. As expected, at the beginning the SNR is 
large and enough to achieve good tracking, but as the time passes the missile moves away 


from the radar, so the SNR decreases. 
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Figure 25. SNR versus Time for the Radar #1 


Figures 26 and 27 represent the SNR performance of the other six radars. Taking 
the position of the sensors relative to the missile into account, the three first sensors 
(RF#1, RF#2, RF#3) have better performance compared to the fourth sensor (RF#4). The 
fifth sensor (RF#5) has the same response as the previous sensors but provides smaller 
values than the first three sensors because it is placed away form the route of the missile. 
For the next two sensors, which are farther away from the launch site, the SNR responses 
are presented in Figure 27. For RF#6 and RF#7, the beginning values of the SNR are very 


small, but as the time passes, the missile approaches the sensors, so the values of the SNR 
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increase. Especially, for the sixth sensor (RF#6), the SNR reaches the peak value at 240 


sec of the flight; then the missile moves away from the sensor, so SNR decreases. 
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Figure 26. 


SNR versus Time from Radars #2, #3, and #4. 
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Figure 27. SNR versus Time from Radars #5, #6, and #7. 


We now present the data fusion problem. The calculation of the track scoring 


function is to identify the sensor with the best track file. A track score is calculated for 


each sensor based on the kinematics of the missile flight parameters and the SNR at the 


sensor. By using likelihood ratios, the optimum track file of the threat can be determined 


and the corresponding track file can be transmitted to the Battle Manager Control in order 


to launch the interceptor vehicle against the threat. Using the optimum track scoring 


techniques developed in this thesis, the best track file can be sent to the interceptor to 


guide it to the ballistic threat. This leads to a faster response and the threat can be 


destroyed inside the territory of the country which launched the missile before the missile 


employs any countermeasures. Figure 28 shows the track score for three RF sensors as a 
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function of time (or scans) and demonstrates that during the threat engagement, some 
sensors have higher track score at certain times while others have a higher track scores at 
other positions of the threat flight. The objective of the track score processing is to 
calculate the track score for each track and select the optimum track that should be used 


to drive the interceptor. 
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Figure 28. Track scores for Radars #1, #2, and #3. 
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Based on the track score in Figure 28 at the beginning of the flight, the third radar 
(RF#3), which is closer to the threat at the beginning, has the best track score quality. As 
the flight continues, the first radar (RF#1) demonstrates the best track score for the next 
fifty seconds. Finally, the second radar (RF#2) has the best track score until the end of the 


boost phase. 
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B. IR/LIDAR PERFORMANCE 


In this section, we present the results of the IR and the LIDAR sensors placed on 


the satellites platforms, SAT#1 and SAT#2. The IR and LIDAR sensors on each platform 


have the same design parameters. In Chapter I, the SNR as a function of time for the IR 


and LIDAR sensors placed on SAT#1 were shown (Figure 9). Figure 29 presents the 
SNR as a function of time for the IR sensors placed on SAT #1 and SAT#2. 
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Figure 29. 


SNR versus time for the IR systems on SAT#1, and SAT#2 


Comparing the two plots in Figure 29 for the two IR sensors on SAT#1 and 


SAT#2, both IR systems show satisfactory performance at the beginning of the flight, 


having approximately the same value of 25 dB. As the time passes, the SNR for both IR 
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sensors increases due to the fact the missile approaches the two platforms. The second 
sensor has better performance as its peak is approximately 42 dB versus 38 dB for the 
first sensor. Then SNR starts to decrease for both sensors as the missile moves away from 


the two platforms. 


The LIDAR sensors are placed on the same platforms as the IR sensors; the SNR 
performance of the LIDAR on SAT#1 was shown in Figure 10. Figure 30 presents the 
SNR of the LIDAR sensors placed on SAT#1 and SAT#2. (The SNR plot of LIDAR #1 


is repeated here from Figure 10). 
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Figure 30. SNR versus time for LIDAR systems on SAT#1, and SAT#2. 


Comparing the two plots in Figure 30, the two systems provide the same 


performance at the beginning of the trajectory as the values of SNR are the same and 
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approximately equal to 20 dB. As the time passes, the second system performs better than 
the first system. The peak of the second system is equal to 51 dB compared 44 dB for the 


second system. 


Also, from Figures 29 and 30, the graphs are similar in shape. Only the values of 
the SNR are different in each graph. Even thought the graphs have been derived from 
different sensors (IR and LIDAR), their tracking responses are the same. Thus, two 


different sensors placed on the same platform confirm the trajectory of the missile. 


In summary, in this chapter, the results from the simulation for the seven ground 
radars as well as the results for the spaceborne IR and LIDAR sensors were presented. 
Also, the results for the data fusion of the radars were presented. Finally, the performance 
of the spaceborne IR and LIDAR sensors is compared for evaluating the accuracy of the 


simulation to detect a ballistic missile in the boost-phase. 
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Vv. CONCLUSION 


In this thesis, a scenario of many different types of sensors in the vicinity of a 
launched ballistic missile is studied. The main goal was the evaluation of the performance 
of these sensors during the boost-phase stage of the missile’s flight. The SNR of the 


sensors is used as a figure of merit. 


Seven X-band, ground based radars with the same parameters are deployed at 
different positions. Observing the results of the radars, it is inferred that the distance 
between the radar and the target is the basic factor to be taken into account for evaluating 


the performance of the sensor in the data fusion process. 


Two IR sensors are deployed in two LEO satellites in order to evaluate their 
performance during boost-phase. They have an advantage over radars due to their 
capability to detect the radiation of the missile in this phase. Two LIDAR sensors are also 
deployed on the same platforms along with the IR sensors. The LIDAR sensor is used to 
provide the third dimension, the range, which an IR sensor cannot calculate. Based on the 
results, it is inferred that the use of the LEO satellite supported by IR sensors improves 
the detection and tracking performance of a ballistic target in the first stage of its flight. 
The drawback is that the area covered by a LEO satellite for surveillance is limited, so a 


network of LEO satellites is required for covering the Earth. 


The next goal of the thesis was the fusion of the radar sensor data. The technique 
of likelihood ratios for tracking of each sensor is used. The evaluation provides results 
about the performance of each sensor with respect to the scan number, which in this work 


represents the time of the flight. 


A. SUGGESTIONS FOR FUTURE WORK 


The implementation of different types of radar in the vicinity of the missile launch 
site can be an issue for further investigation. Radars with different parameters, such as 


transmitted power, frequency, PRF, antenna gain and antenna temperature based on 
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actual sky temperatures for the radar frequency and beam elevation angles of interest will 
present a more realistic scenario for the data fusion problem. In this situation, the distance 
between the sensor and the target is not the primary factor for evaluating the data fusion 


problem. 


The data fusion problem for the IR sensors based on the likelihood ratio can be 
studied in a future research. The fusion of the IR sensor data in boost phase is not the 
primary goal as the function of IR is only to detect the ballistic missile. In the next stages 
is where the IR sensor can contribute to the tracking of the missile hence the data fusion 


problem. 


A study involving deployment of multiple simultaneously launched ballistic 
missiles is of interest. A network of LEO satellites with IR sensors may be employed to 


discriminate targets and determine the tracks of multiple targets. 
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APPENDIX. FLOW CHART CODE 


The MATLAB code used in this work is based on the code of Reference [2]. The 
code has been modified in order to include the function of the data fusion among the 


radars and the deployment of the IR and LIDAR sensors. Figures 31, 32 and 33 present 


the flow charts for data fusion. 
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Figure 31. RFsimulation() for the seven radars 
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Figures 34 and 35 present the flow chart for the LIDAR sensors 
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Figure 34. LIDARsimulation() for the two LIDAR sensors 
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Figure 36 presents the flow chart for IR sensors. 
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Figure 36. Flow chart for IR sensors 
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